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ABSTRACT 

We present the first absolute proper motion measurement of Leo I, based on two epochs of Hubble 
Space Telescope (HST) ACS/WFC images separated by ~ 5 years in time. The average shift of Leo I 
stars with respect to ~ 100 background galaxies implies a proper motion of (hw,^n) — (0.1140 ± 
0.0295, —0.1256 ± 0.0293)mas yr~ 4 . The implied Galactocentric velocity vector, corrected for the 
reflex motion of the Sun, has radial and tangential components Vrad = 167.9 ± 2.8 kms -1 and Vt a n = 
101.0±34.4kms~ 1 , respectively. We study the detailed orbital history of Leo I by solving its equations 
of motion backward in time for a range of plausible mass models for the Milky Way and its surrounding 
galaxies. Leo I entered the Milky Way virial radius 2.33 ± 0.21 Gyr ago, most likely on its first infall. 
It had a pericentric approach 1.05 ± 0.09 Gyr ago at a Galactocentric distance of 91 ± 36 kpc. We 
associate these time scales with characteristic time scales in Leo I's star formation history, which shows 
an enhanced star formation activity ~ 2 Gyr ago and quenching ~ 1 Gyr ago. There is no indication 
from our calculations that other galaxies have significantly influenced Leo I's orbit, although there is 
a small probability that it may have interacted with either Ursa Minor or Leo II within the last ~ 1 
Gyr. For most plausible Milky Way masses, the observed velocity implies that Leo I is bound to the 
Milky Way. However, it may not be appropriate to include it in models of the Milky Way satellite 
population that assume dynamical equilibrium, given its recent infall. Solution of the complete (non- 



radial) timing equations for the Leo I orbit implies a Milky Way mass M\ 



MW ,vir 



3.15 



+ 1.58 
-1.36 
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with the large uncertainty dominated by cosmic scatter. In a companion paper, we compare the new 
observations to the properties of Leo I subhalo analogs extracted from cosmological simulations. 
Keywords: Astrometry — galaxies: kinematics and dynamics — galaxies: individual (Leo I) — Galaxy: 
kinematics and dynamics — Galaxy: halo — Local Group 



1. INTRODUCTION 

Structures in the Universe cluster on various scales. 
The Milky Way (MW) is no exception, as evidenced by 
its system of satellites. Most of these satellites are dwarf 
spheroidal galaxies (dSphs). These are the most dark 
matter dominated stellar systems currently known, with 
mass-to-light (M/L) ratios of up to a few thousand in 
units of M /L Q (e.g., Wolf et al. 2010) 

In the current paradigm for galaxy formation, dark ha- 
los of galaxies form through the accumulation of smaller 
subunits. The MW satellite system is one of the best ob- 
jects to study these hierarchical evolutionary processes 
in action, due to its proximity. In the last decade, many 
wide-field ground-based surveys have led to discoveries 
in this area. For example, the The Two Micron All Sky 
Survey (2MASS) unveiled the ongoing disruption of the 
Sagittarius dSph that has produced a giant stream of 
stars wrapping around the entire MW at least a single 
time (e.g., Majewski et al. 2003), and the SDSS has re- 
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vealed many other such streams that once belonged to ei- 
ther dwarf galaxies or globular clusters (Grillmair 2009, 
and references therein). In addition, there is evidence 
for recent accretion and build up of the MW satellite 
system: HST proper motion measurements of the two 
most-massive MW satellites, the Large and Small Magel- 
lanic Clouds (LMC and SMC) (Kallivayalil et al. 2006a,b; 
Piatck et al. 2008), suggest that these galaxies were not 
born as MW satellites but instead may be falling into 
the Local Group for the first time (Besla et al. 2007; 
Kallivayalil et al. 2012). 

As tracer objects, MW satellites are valuable tools for 
studying the size and mass of the MW halo because their 
orbits contain important information about the host po- 
tential. Distant satellite galaxies are of particular inter- 
est because: (1) they probe the dark halo at the largest 
radii; and (2) their kinematics may not have been fully 
virialized yet. Measuring the space motions of distant 
satellites with respect to the MW is therefore crucial for 
gaining insights into the MW virial mass and the mass 
assembly at late epochs. 

So far, there are only three known objects thought to 
be associated with the MW at a distance beyond 200 
kpc from the Galactic center: the dSphs Leo I, Leo II, 
and Canes Venatici I. Leo I, unlike the others, has an 
unusually large Galactocentric radial velocity at its ex- 
treme distance (Mateo et al. 1998, 2008). Because of 
this, Leo I has played an important role in our interpre- 
tation of the MW satellite system. One reason for this is 
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that Leo I disproportionately affects MW mass estimates 
based on the assumption of equilibrium kinematics: in- 
cluding or excluding it from the MW satellite population 
sample produces very different estimates (e.g., Zaritsky 
et al. 1989; Kulessa & Lynden-Bell 1992; Kochanek 1996; 
Wilkinson & Evans 1999; Watkins et al. 2010). 

Another issue on which Leo I has generated much inter- 
est and debate is on the topic of its specific orbit. Its large 
radial velocity led Byrd et al. (1994) to suggest that Leo I 
once belonged to M31 and is now on a hyperbolic orbit 
flying past the MW. However, in an earlier study, Zarit- 
sky et al. (1989) argued against the possibility that Leo I 
is unbound to the MW. More recently, Sohn et al. (2007) 
and Mateo et al. (2008) carried out orbital analyses com- 
bined with high-precision radial velocities of individual 
Leo I member stars to study the orbit in detail. The 
former study suggested that Leo I was tidally disrupted 
on one or two perigalactic passages about a massive Lo- 
cal Group member (most likely the MW), whereas the 
latter study proposed involvement of a third body that 
may have injected Leo I into its present orbit a few Gyr 
before its last perigalacticon. The orbit of Leo I can also 
shed light on studies using satellites as test particles in 
a cosmological context (e.g., Li & White 2008; Rocha et 
al. 2012). Unfortunately, many orbital scenarios remain 
possible as long as only one component of Leo I's velocity 
(along the line of sight) is known. To make progress, it 
is necessary to know also the proper motion of Leo I to 
yield the full three-dimensional Galactocentric velocity. 

Due to its distance, it has not previously been possible 
to measure the proper motion of Leo I. The most distant 
MW satellite with a measured proper motion so far is 
Leo II, for which a measurement was obtained with the 
second-generation HST instrument WFPC2 (Lepine et 
al. 2011). However, the large uncertainty of this proper 
motion measurement with an accuracy of 0.132masyr~ 1 , 
corresponding to 144 km s -1 at the distance of Leo II, 
limits its usefulness in constraining models. Leo I is lo- 
cated at a distance of <~ 260 kpc, which is even <~ 40 kpc 
farther away than Leo II. 

However, we have recently pioneered a method to 
measure the proper motions of galaxies as far away as 
M31 (Sohn et al. 2012) using the third- and/or fourth- 
generation HST instruments ACS and WFC3. This in- 
volves sophisticated data analysis techniques to measure 
from deep images taken years apart the relative shifts 
of thousands of stars in the galaxy with respect to hun- 
dreds of distant compact galaxies in the background. We 
were able to achieve a final proper motion accuracy of 
0.012 masyr -1 , which yields a velocity uncertainty of 
44kms _1 at the distance of M31, using 18 independent 
measurements on three different fields. Leo I is at a dis- 
tance three times closer than M31, so the same level of 
velocity uncertainty is well within the reach of our proven 
techniques. Because this would certainly yield useful new 
constraints on the orbit of Leo I, we designed an obser- 
vational program to measure the absolute proper motion 
of Leo I for the first time. We report here on the results 
and implications of our program. 

This paper is organized as follows. In Section 2 we 
describe the data and analysis steps, and we present the 
inferred Leo I proper motion. In Section 3, we correct 
the measured proper motion for the solar motion to de- 



rive the Galactocentric space motion. We then explore 
the implications for the past orbit of Leo I under a vari- 
ety of assumptions for the mass and mass distribution of 
the MW and other Local Group galaxies. In Section 4, 
we discuss the implications of the results for our under- 
standing of both Leo I and the MW, and summarize the 
main results of the paper. 

This is the first paper in a series of two. In Paper II 
(Boylan-Kolchin et al. 2012b), we compare the new ob- 
servations of Leo I to the properties of Leo I subhalo 
analogs extracted from state-of-the-art cosmological sim- 
ulations. We use this comparison to place additional 
constraints on the mass of the MW, the properties of 
its satellite system, and the past history of Leo I. 

2. THE PROPER MOTION OF LEO I 
2.1. Hubble Space Telescope Data 

The data used in this study to measure the proper 
motion of Leo I consist of images taken with HST in two 
different epochs separated by ~ 5 years in time. For the 
first-epoch data, we used the images taken in February 
2006 for the science program GO-10520 (PI: T. Smecker- 
Hane) to study the star formation history of Leo I. A 
field slightly offset from the center of Leo I was imaged 
with the ACS/WFC using the F435W and F814W filters. 
Figure 1 shows the field location overlaid on a STScI 
Digital Sky Survey image centered on Leo I. Sets of 7 
and 6 images with exposure times of 1,700 sec each were 
obtained in F435W and F814W, respectively, and three 
additional 440 sec F814W images were obtained as well. 

The second-epoch data were obtained in January 2011 
for our science program GO-12270 (PI: S. Sohn). We pre- 
analyzed the first-epoch data to enable optimal design of 
the second-epoch observations. This analysis indicated 
that the F814W filter provides a slightly better astro- 
metric handle on extended objects than the F435W filter 
(see also Mahmud & Anderson 2008, for a discussion of 
the wavelength dependence of the astrometric accuracy) . 
We therefore took second-epoch observations only with 
F814W. We obtained 12 images with individual exposure 
times ranging from 1,267 to 1,364 sec. The resulting total 
exposure time for the second epoch was <~ 16 ksec. Indi- 
vidual exposures were dithered using a pattern designed 
to optimize the sampling of the point-spread functions 
(PSFs) for stars that fall on different parts of the detec- 
tor. This "pixel-phase coverage" is crucial for creating a 
high-resolution stacked image from a limited number of 
individual exposures. 

We matched the orientation and field center of the 
second-epoch observations as closely as possible to those 
of the first-epoch observations. However, due to unavail- 
ability of the same guide stars used for the first-epoch 
observations, we had to use an orientation that differed 
by ~ 1°. We also obtained parallel observations with 
WFC3/UVIS in the second epoch for an off-target field. 
This field, also shown in Figure 1, was imaged in F438W 
and F814W to allow a study of the stellar population in 
the outer halo of Leo I. However, these observations are 
not discussed further in the present paper. 

2.2. Measurement Technique 

To measure the proper motion of Leo I we compare 
the two epochs of HST imaging data, and determine the 
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Figure 1. Field location of our ACS/WFC observations plotted in red over a 15'xl5' section of the sky centered on Leo I from the STScI 
Digital Sky Survey. The dashed ellipse represents the King limiting boundary as derived by Sohn et al. (2007). The parallel WFC3/UVIS 
field is plotted in blue dashed lines. 



average shift of Leo I stars relative to distant background 
galaxies. This requires a method that allows accurate po- 
sitions to be measured for both stars and compact galax- 
ies. Mahmud & Anderson (2008) presented a method 
that accomplishes this by constructing and fitting an in- 
dividual template for each source in an image. Sohn et al. 
(2012) implemented, expanded and applied this method 
to measure the proper motion of the galaxy M31 using 
HST ACS/WFC and WFC3/UVIS images of three dis- 
tinct fields imaged over a 5-7 year time baseline. We 
adopt their method here to also analyze the new Leo I 
data. We discuss only the main outline and results of the 
proper motion derivation, and refer the reader to Sohn et 
al. (2010) and Sohn et al. (2012) for more details about 
the methodology. 

All the science _f It .fits images for the first and sec- 
ond epochs were downloaded from the archive. To each 
image we first applied the CTE-correction routine devel- 
oped by Anderson & Bedin (2010). We then used the 
img2xym_WFC . 09x10 program (Anderson & King 2006) 
to determine a position and a flux for each star in each 
exposure. The positions were subsequently corrected for 
the known ACS/WFC geometric distortions. Separate 
distortion solutions were used for the first- and second- 
epoch data, to account for a difference between pre- and 
post-SM4 ACS/WFC data (see Section 3.3 of Sohn et al. 
2012). We then adopt the first exposure of the second 
epoch (jbjmOlkkq) as the frame of reference. We cross- 
identify stars in this exposure and the same stars in other 
exposures. We use the distortion-corrected positions of 
the cross-identified stars to construct a six-parameter lin- 
ear transformation between the two frames. These trans- 
formations are then used in a program that constructs a 
stacked image, cleaned of cosmic rays and detector ar- 
tifacts, of the different exposures in each filter+epoch 



combination. The stacked images were super-sampled 
by a factor of 2 relative to the native ACS/WFC pixel 
scale for better sampling. 

Stars and galaxies were identified from the stacked 
second-epoch F814W image, which provides our deep- 
est view of the field. First, a list of point 
sources was constructed from the sources detected by 
img2xym_WFC . 09x10. The selection of Leo I stars from 
this list is relatively straightforward, since our target field 
is numerically and spatially dominated by Leo I stars 
(Leo I is located at high Galactic latitude b = 49°, so 
few foreground stars are expected). To select only bona 
fide and well-measured Leo I stars, we require: (1) small 
RMS scatter between the 12 independent position mea- 
surements; (2) consistent position in the color-magnitude 
diagram (obtained from combination with the first-epoch 
F435W results) with the expected Leo I stellar evolu- 
tionary features; and (3) consistent proper motion with 
the other Leo I stars. This yielded a list of 36,000 stars 
suitable for proper motion analysis. For the selection of 
background galaxies in the target field, we started with a 
candidate list generated by running SExtractor (Bertin 
& Arnouts 1996) on the stacked image. From this list we 
then carefully identified 116 compact background galax- 
ies by eye. 

For each star/galaxy in each F814W exposure of each 
epoch we then measured in consistent fashion a posi- 
tion using the template-fitting method. Templates were 
constructed from the second-epoch stacked image via in- 
terpolation. The template-fitting for the first-epoch data 
included an additional 7 x 7-pixel convolution kernel, to 
allow for PSF differences between epochs. This kernel 
was determined from the data for bright and isolated 
Leo I stars, without any assumed field-dependence. This 
is similar to what was done in our analysis of M31 (Sohn 
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Figure 2. Onc-dimcnsional positional uncertainty per exposure as 
a function of instrumental magnitude, for Leo I stars (upper panels) 
and background galaxies (lower panels). The error is defined as 

(Tl_£> = y \(°~x 2 + o-y 2 ). Here a x and a y arc the pcr-coordinatc 

RMS residuals with respect to the average for the multiple second- 
epoch measurements. 

et al. 2012), except that in that case the role of the epochs 
was reversed (since for M31 our first-epoch data were the 
deepest). 

The template-fitted positions were corrected as before 
for the ACS/WFC geometric distortion. Again, star po- 
sitions in individual exposures were used to determine 
six-parameter linear transformations with respect to the 
first exposure (jbjmOlkkq, adopted as the frame of refer- 
ence) but now using only the selected Leo I stars. These 
linear transformations were then used to transform the 
measured positions of all selected stars and background 
galaxies in all exposures into the reference frame. 

The individual exposures lead to multiple determina- 
tions for the position of each star or background galaxy 
in each epoch. We average these determinations to ob- 
tain the average position of each source in each epoch. 
The RMS scatter between measurements quantifies the 
random positional uncertainty in a single measurement. 
The error in the mean is smaller by y/~N, where N is the 
number of exposures in the epoch. 

Figure 2 shows the one-dimensional positional errors 
per second-epoch exposure for Leo I stars (upper pan- 
els) and background galaxies (lower panels), as a func- 
tion of instrumental magnitude (defined as mi nstr = 
— 2.5 log[counts]). Brighter objects produce a higher 
signal-to-noise ratio than fainter objects, and therefore 
have more accurately determined positions. Stars are 
more compact than background galaxies, and therefore 
generally have more accurately determined positions. 
However, the brightest galaxies have positional uncer- 
tainties comparable to those of the brightest stars. 

The proper motion of a source is the difference between 



its average position in the two epochs, divided by the 
time baseline (4.93 years). By construction, our method 
aligns the star fields between the two epochs. So with 
this convention, Leo I stars have zero motion on average, 
and the background galaxies move over time. Of course, 
in reality the background galaxies are stationary and the 
Leo I stars move in the foreground. So the actual average 
proper motion of Leo I stars is obtained as —1 times 
the measured average proper motion of the background 
galaxies. 

Figures 3 shows the measured proper motion of each 
star in X and Y as a function of detector coordinates, 
for one of the 1,700 sec first-epoch images. The motions 
are zero on average by construction, but there remain 
small residual trends with position on the detector at lev- 
els < 0.01 pixel. This could be due, e.g., to limitations 
in the adopted geometric distortion corrections. These 
trends are corrected by measuring the displacement of 
each background galaxy with respect to only those Leo I 
stars that lie in the vicinity of the galaxy. This "lo- 
cal correction" removes any remaining systematic proper 
motion residuals associated with the detector position. 
Each local correction was constructed using stars of sim- 
ilar brightness (±1 mag) and within a 200 pixel region 
centered on the given background galaxy. 

Figure 4 shows the measured proper motion of each 
background galaxy in X and Y as a function of detector 
coordinates, for the same 1,700 sec first-epoch image as in 
Figure 3. The proper motion of each background galaxy 
is measured with respect to the average second-epoch 
position of that galaxy, including the local correction. 
There are far more stars than background galaxies in our 
images, and star positions are generally determined more 
accurately than galaxy positions. Hence, the final proper 
motion uncertainty is dominated by the astrometric ac- 
curacy for the background galaxies. For each individual 
first-epoch exposure we take the weighted average over all 
background galaxies (larger symbols in Figure 4 denote 
background galaxies with more accurate measurements, 
which receive more weight) to obtain a individual Leo I 
proper motion estimate (red line). The 1-cr confidence 
region around the weighted average (dashed lines) was 
computed using the bootstrap method (Efron & Tibshi- 
rani 1993), with 10,000 bootstrap samples. 

2.3. Inferred Proper Motion 

The proper motion diagram for the 9 independent first- 
epoch measurements is shown in Figure 5. We plot the 
data points for the longer (1,700 sec; open squares) and 
shorter (440 sec; open triangles) exposures in the same 
diagram. We transformed the proper motions and their 
associated errors along the detector axes to the directions 
West and North using the orientation of the reference 
image with respect to the sky (— 47?7). Table 1 lists the 
proper motion for each first-epoch image and the cor- 
responding error, along with the number of background 
galaxies used for the proper motion derivation. 

The final average proper motion of Leo I is calculated 
by taking the error-weighted mean of the 9 independent 
measurements listed in Table 1. This yields 

(Vw^n) = (0.1140±0.0295,-0.1256±0.0293)mas yr" 1 . 

' (1) 

This result differs from zero at approximately 4-cr confi- 
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X(pbe) Y(plx) 

Figure 3. Displacements of individual stars (dark gray dots) versus detector location between one of the 1,700 sec first-epoch images 
(j9gz04tsq) and the average of the second-epoch images, plotted separately for x and y positions. The average and RMS displacements 
of stars in every 500-pixel bin are shown in red. The displacements average to zero by construction. The scatter is a measure of the 
per-exposure positional accuracy for a star of average brightness. Low-level trends are indicative of residual detector effects. The units arc 
in native ACS/WFC pixels, and X and Y positions are in the reference frame. 
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Figure 4. Displacements of background galaxies versus detector location between one of the 1,700 sec first-epoch images (j9gz04tsq) 
and the average of the second-epoch images, plotted separately for X and Y positions. The black points show the relative displacements 
measured for different background galaxies. The weighted average for all galaxies is shown as the red line; dashed red lines indicate the l-a 
error region around the average. This region is smaller than the scatter between the points by a factor of ~ V^V, where N is the number 
of background galaxies. The radius of each black point is proportional to 1/A , where A is the proper motion measurement uncertainty 
for the particular background galaxy. Hence, the area of each point is proportional to the weight a point receives in the final weighted 
average. Symbols in green in the top left panel illustrate how the point size relates to the proper motion uncertainty A. The units are in 
native ACS/WFC pixels, and X and Y positions are in the reference frame. 
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Figure 5. Leo I proper motion results. Each gray dot with an 
error bar indicates the average proper motion of Leo I stars inferred 
from a single first-epoch exposure. Measurements using images 
with longer (closed square) and shorter (closed triangle) exposure 
times are indicated with different symbols. The solid red data point 
is the weighted average of the nine separate measurements, which 
is the final result of our analysis. The origin corresponds to the 
velocity such that Leo 1 has no transverse motion in the heliocentric 
rest frame. The star symbol corresponds to the velocity such that 
Leo I has no tangential velocity in the Galactocentric rest frame 
(i.e., a radial orbit with respect to the MW). 



dence in each coordinate direction, so the detected mo- 
tion of Leo I is very statistically significant. 
The quantity 



X 



E 



( Be. 



(2) 



provides a measure of the extent to which different mea- 
surements agree to within the random errors. In ab- 
sence of systematic errors, one expects that this quan- 
tity follows a x 2 probability distribution with Njjp = 
18 — 2 — 16 degrees of freedom. The expectation value 
for such a distribution is NdFi an d the dispersion is 
~ \/2Ndf = 5.7. We find \ 2 = H-2 for our measure- 
ments. This indicates that the measurements from the 
different exposures are consistent, and that the errors 
may actually be slight overestimates. 

The final Leo I proper motion uncertainties correspond 
to ~ 29 fj,a,s yr" 1 . This is a factor of > 2 larger than what 
we achieved for M31 (Sohn et al. 2012), due mostly to 
the fact that for M31 deeper exposures were available 
for three different fields. Our M31 HST measurements 
approached the accuracy achieved using VLBA water 
maser observations for the M31 satellites M33 and IC10. 
Whereas this is not the case here, our Leo I measure- 
ments are more accurate than what has been achieved 
with HST for other MW satellites, using one or more 
fields centered on background quasars (see the compila- 
tion in Table 4 of Watkins et al. 2010). 

3. THE ORBIT OF LEO I 
3.1. Velocity in the Galactocentric Rest Frame 



Table 1 

Proper Motion Results 
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3gz05ulq 


0.0528 ± 0.0996 


-0.0865 ± 0.0962 
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-0.0966 ± 0.0970 
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3gz06ksq 


0.0545 ± 0.0917 


-0.0718 ± 0.0939 


53 
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3gz06kuq 


0.0501 ± 0.0894 


-0.1449 ± 0.0891 


54 


w 


eighted av. b 


0.1140 ± 0.0295 


-0.1256 ± 0.0293 





a Number of background galaxies used for deriving the average 
proper motion of each field. 

b Weighted average of the results for the 9 independent measure- 
ments. 



We adopt a Cartesian Galactocentric coordinate sys- 
tem (X,Y,Z), with the origin at the Galactic Center, 
the X-axis pointing in the direction from the Sun to the 
Galactic Center, the V-axis pointing in the direction of 
the Sun's Galactic rotation, and the Z-axis pointing to- 
wards the Galactic North Pole. The position and velocity 
of an object in this frame can be determined from the ob- 
served sky position, distance, line-of-sight velocity, and 
proper motion, as in, e.g., van der Marel et al. (2002). 

To determine the Galactocentric position of an object, 
it is necessary to also know the distance Rq of the Sun 
from the Galactic Center. Moreover, it is necessary to 
know the velocity of the Sun inside the MW to turn 
observed heliocentric rest-frame velocities into Galacto- 
centric rest-frame velocities. Following van der Marel 
et al. (2012a), we adopt the recent values of McMillan 
(2011) for the distance of the Sun from the Galactic 
center and the circular velocity of the local standard of 
rest (LSR): R = 8.29 ± 0.16 kpc and V = 239 ± 5 km 
s _1 . For the solar peculiar velocity with respect to the 
LSR we adopt the estimates of Schonrich et al. (2010): 
(Ppec, V poc W pcc ) = (11.10, 12.24, 7.25) km s" 1 with un- 
certainties of (1.23,2.05,0.62) km s" 1 . 

To obtain the distance of Leo I, we average the dis- 
tances measured via the tip of the red-giant branch 
(TRGB) method in the last decade (Mendez et al. 2002; 
Bellazzini et al. 2004; Held et al. 2010), which yields 
256.7 ± 13.3 kpc. This implies a Galactocentric (X, Y, Z) 
position 



-125.0, -120.8, 194.1) kpc, 



(3) 



with an uncertainty of 13.3 kpc along the line-of-sight 
direction. 

The most recent measurement of the systemic heliocen- 
tric line-of-sight velocity of Leo I is t>Los — 282.9 ± 0.5 
km s _1 (Mateo et al. 2008). The measured proper mo- 
tion \x from equation (1) corresponds to a heliocentric 
transverse velocity in km/s equal to 4.7404 x £>(kpc) x 
/i(mas yr _1 ). This implies 

(v w ,v N ) = (138.7 ±36.6, -152.8 ± 36.5) kms" 1 (4) 

with proper motion errors dominating over distance er- 
rors in the determination of the velocity errors. The in- 
ternal velocity dispersion of Leo I is u = 9.2±0.4kms~ 1 , 
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with little evidence for rotation (Mateo et al. 2008). This 
is well below our observational velocity errors. Hence, 
there is no need to correct the observed values for the 
internal kinematics of Leo I, even though our field was 
offset from its photometric center (see Figure f). 

The velocity for which Leo I would be on a radial or- 
bit with respect to the MW (i.e., the velocity for which 
there is zero tangential velocity in the Galactocentric rest 
frame) is 

(vw,v N ) iad = (58.0, -214.0) kms" 1 . (5) 

This differs from the measured proper motion at almost 
3-ct significance. Therefore, our measurements imply 
that Leo I is not on a radial orbit about the MW. 

Several authors have argued previously that the Galac- 
tocentric tangential velocity of Leo I is probably small, 
given its significant radial velocity (Byrd et al. f994; 
Sohn et al. 2007; Mateo et al. 2008). Indeed, Figure 5 
shows that the observed proper motion does fall in the 
same quadrant of proper motion space as a radial or- 
bit. This can be interpreted as a consistency /plausibility 
check on the proper motion measurement. The same was 
found for the case of M31, for which we also presented 
several other successful consistency checks on our proper 
motion analysis methodology (Sohn et al. 2012). 

Conversion into the Galactocentric rest frame yields 
for the velocity vector of Leo I 

v= (-167.7±31.9, -37.0±33.0,94.4±24.2)kms _1 . (6) 

The listed uncertainties here and hereafter were obtained 
from a Monte-Carlo scheme that propagates all observa- 
tional distance and velocity uncertainties and their cor- 
relations, including those for the Sun. Note that the 
Galactocentric velocity uncertainties are highly corre- 
lated because the observational velocity uncertainty is 
much larger in the transverse direction than in the line- 
of-sight direction. 

The corresponding Galactocentric radial and tangen- 
tial velocities are 

(Viad.Vtan) - (167.9 ± 2.8, 101 .0 ± 34.4) km s- 1 . (7) 

Although the tangential velocity is significantly non-zero, 
it is less than the radial velocity. So while Leo I is not 
on a radial orbit about the MW, the orbit must be fairly 
elliptical. The observed total Leo I velocity with respect 
to the MW is 

v = \v\ = 196.0 ± 19.4 km s" 1 (8) 

with the listed numbers corresponding to the peak and 
the symmetrized lc of the v probability distribution. 6 

These inferred Leo I velocities use a solar velocity in- 
side the MW based on Schonrich et al. (2010) and McMil- 
lan (2011), which yields an azimuthal velocity compo- 
nent i>0,© = 251.2 kms . However, alternative values 
for the solar velocity continue to be in common use. 
These differ from the values used here primarily in the 

6 The error distribution of the total velocity v is somewhat asym- 
metric, but this is more pronounced at the 2cr than at the la level. 
The median of the v distribution and the surrounding 68% (95%) 

confidence intervals are v = 199.8^^ 1(^293) kms -1 , as used in 
Paper II. 



azimuthal direction. For example, with the old IAU rec- 
ommended circular velocity Vo = 220 kms -1 and the pe- 
culiar velocities from Dehnen & Binney (1998), v^q = 
225.2 kms -1 . Based on these latter values, the Galac- 
tocentric Leo I velocities would be V Ta d = 180.5 kms -1 , 
Vtun = 93.8 kms -1 , and v = 203.4 kms -1 , which can 
be compared to the values in equations (7) and (8). 
While the change in T^ad is significant compared to the 
small uncertainties, V tan and v change by much less than 
the observational uncertainties. The conclusions of the 
present paper and Paper II are therefore not very sen- 
sitive to the adopted solar velocity. The value of v^.© 
assumed here is consistent with the recent determination 
of W0, = 242± 10 kms" 1 by Bovy et al. (2012), while the 
alternative W0 i0 = 225.2 kms -1 discussed above is not. 7 

3.2. Keplerian Orbit Calculations 

To assess the implications of the new measurements, 
we start with the assumption that the MW can be ap- 
proximated as a point mass, and that Leo I orbits in its 
potential as a test particle on a Keplerian orbit. The as- 
sumption of a Keplerian potential for the MW is not as 
unreasonable as it may seem at first. The large Galacto- 
centric distance of Leo I, r = |f| = 260.6 kpc, combined 
with its significant tangential velocity, implies that much 
of the MW's mass is inside the Leo I orbit at all times. 
We calculate models with more realistic MW potentials 
in Section 3.3. 

The escape velocity for a point mass Mmw is 

v e = V2GM M wA- (9) 

Hence, given the observed total Leo I velocity with re- 
spect to the MW given by equation (8), Leo I is bound to 
the MW if M MW > (1.16 ± 0.24) x 10 12 M Q . Cosmolog- 
ical simulations imply that it is very unlikely to find an 
unbound satellite at the present epoch near a MW-type 
galaxy (e.g., Deason et al. 2011, and Paper II). So if we 
assume that Leo I must be bound to the MW, then this 
can be interpreted as a new crude lower limit on the MW 
mass. 

Alternatively, one may assume that the mass of the 
MW is already constrained from other arguments. In 
that case, one can use the new measurement to assess 
the probability that Leo I is in fact bound. Studies 
of the MW mass have advocated many different val- 
ues, roughly covering the range 0.75-2.25 x 10 12 M Q (see 
Boylan-Kolchin et al. 2012a, for a compilation of recent 
mass estimates of the MW). van der Marel et al. (2012a) 
assumed a flat prior probability over this range, and then 
used a Bayesian scheme to include the latest measure- 
ments of the Local Group timing mass, based on our 
M31 HST proper motion work. The Local Group tim- 
ing mass is relatively high, which increases the likelihood 
of high MW masses compared to low MW masses. We 
combined the probability distribution for the MW mass 
from Figure 4 of van der Marel et al. (2012a) with the 
measured v for Leo I. This implies that there is a 77% 
probability that Leo I is bound to the MW, and a 23% 
probability that it is not bound. The preference for a 

7 Bovy et al. (2012) advocate for a lower circular velocity than 
McMillan (2011), but the only quantity that matters for the cal- 
culations presented here is «<£,©. 
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bound state is consistent with expectations from cosmo- 
logical simulations (e.g., Benson 2005; Wetzel 2011). 

For any assumed point mass Mmw, and given Galac- 
tocentric Leo I phase-space vectors r and v, the shape 
of the Keplerian orbit is determined analytically. We 
calculated these orbits in a Monte-Carlo sense. At each 
Montc-Carlo step we draw a mass Mmw from the pre- 
viously discussed probability distribution derived by van 
der Marel et al. (2012a), and we draw Leo I phase-space 
vectors from the observationally determined values and 
uncertainties. We then determine the statistics of the 
orbital characteristics over the Monte-Carlo sample. 

The Monte-Carlo analysis yields an average ratio of 
pericenter to apocenter distance r por i/r apo = 0.06 ± 0.03 
for the bound elliptical Keplerian orbits (i.e., orbits with 
eccentricity of less than 1). Leo I has a positive radial 
velocity, and is therefore past pericenter. The pericen- 
tric passage occurred at i por i = 1.18 ± 0.45 Gyr ago at 
a Galactocentric distance of r pcr i = 67 ± 39 kpc. The 
velocity at pericenter was w por i = 561 ± 475 kms -1 . The 
uncertainties in these quantities are determined largely 
by the uncertainties in the Leo I phase-space vectors, 
and much less so by uncertainties in MW mass. For 
example, if Mmw is kept fixed at 1.5 x 10 12 M Q for all 
Monte Carlo drawings, then the orbital characteristics 
become: r pori /r apo = 0.04 ± 0.02, t pori = 1.08 ± 0.13 
Gyr, r pori = 69 ± 38 kpc, w pori = 518 ± 417 kms" 1 . In 
Section 3.3 we compare these simple Keplerian results 
to the results from orbit calculations in more detailed 
cosmologically-motivated halo models. 

Another useful application of Keplerian orbits is 
through the timing argument (Kahn & Woltjer 1959). 
This argument assumes that bound galaxy pairs follow a 
Newtonian Keplerian trajectory starting soon after the 
Big Bang, which corresponds to the "first pericenter." 
The galaxies initially move away from each other due to 
the expansion of the Universe, but then fall back towards 
each other due to gravity. In this picture, Leo I is just 
passed its second pericenter. In general, there are four 
observables (the time since Big Bang i, relative distance 
r, radial velocity V^d, and tangential velocity Vtan) and 
four independent orbital parameters (eccentric anomaly 
rj, semi-major axis length a, eccentricity e, and the to- 
tal mass M). Hence, the Keplerian orbit can be solved 
for analytically, as described in e.g., van der Marel &; 
Guhathakurta (2008). One may call this the "complete 
timing argument" (eta). In many applications however, 
the transverse velocity V tan is not known and it is then 
often assumed that V tan = and e = 1. This yields the 
so-called "radial-orbit timing argument" (rta). 

The timing argument has traditionally been applied to 
the MW-M31 system (see van der Marel et al. 2012a, and 
references therein) , but it can also be applied to the MW- 
Lco I system (Zaritsky et al. 1989). The radial-orbit 
timing argument as applied to Leo I with the previously 
derived Galactocentric position and velocity implies a 
mass M M w,rta = (1.50±0.12) x 1O 12 M (consistent with 
the value previously inferred by Li & White 2008, using a 
slightly different assumed solar velocity) . Any tangential 
velocity increases the timing mass. Since we have now 
measured the tangential velocity of Leo I, we can use 
instead the complete timing argument without assuming 
a radial orbit. This implies a mass Mmw, eta = (1-93 ± 



0.42) x 10 12 M Q . 

The error bars in the listed timing masses reflects only 
the propagation of errors in the observational quantities. 
However, it is important to also quantify any inherent bi- 
ases and cosmic scatter, and to calibrate the timing mass 
to more traditional measures of mass, such as the virial 
mass. Li & White (2008) addressed these issues for the 
radial-orbit timing argument using the cosmological Mil- 
lennium simulation. They identified a set of host galaxy 
- satellite pairs with properties similar to the MW-Leo I 
pair. For these pairs they studied the statistics of the 
ratio M2oo/M rta - In the following we find it more conve- 
nient to use the virial mass M v ; r rather than the quantity 
M200, so where necessary we transform the latter to the 
former using M v i r /M 2 oo = 1-19 (as appropriate for an 
NFW halo of concentration c = 9.5; see the Appendix 
of van der Marel et al. 2012a, for the relevant equations 
and mass definitions). 

Yang-Shyang Li kindly made available the catalog of 
galaxy pairs used in the analysis of Li & White (2008) 
(the sample defined by the bottom row of their Table 3). 
This allowed us to perform an analysis similar to theirs, 
but now for the complete timing argument. We find that 
the bias M v - lr /M ct& = lAQ^^- This estimate was ob- 
tained, as in the Li & White (2008) analysis, by averaging 
over all satellites in the simulation sample, independent 
of tangential velocity. However, we have now measured 
the tangential velocity of Leo I. So we can get a more ap- 
propriate measure of the bias by including only the satel- 
lites with tangential velocities similar to that of Leo I. If 
we require agreement in Vtan to within 25kms _1 , we find 
that M vir /M cta = 1.63t°;^. For comparison, this same 
selection yields for the radial-orbit timing argument that 
M v i r /M rta = 2.10^g'gg. So the complete timing argu- 
ment yields estimates of M v j r that are biased low, but 
not by as much as the radial orbit timing argument. This 
is because part of the bias is due to the fact that satel- 
lite galaxies generally have non-zero tangential velocities, 
and this is explicitly taken into account in the complete 
timing argument. 

These results for cosmic bias and scatter can be com- 
bined with the previously inferred values of Mmw, eta and 
M MW ,rta. This yields M MW ,vir = 3.15±1;|| x 10 12 M Q 
from the complete timing argument, and MMw,vir = 
3.14^1.06 x 10 12 M Q from the radial-orbit timing argu- 
ment, respectively. So the two timing arguments give 
similar results and uncertainties. The mass estimates 
are higher than most MW mass estimates based on other 
methods (consistent with the results of Li & White 2008), 
but they are probably not inconsistent with other MW 
mass estimates given the significant cosmic scatter. This 
situation is similar to what was found for MW mass esti- 
mates based on the timing of the MW-M31 system (van 
der Marel et al. 2012a). 

3.3. Detailed Orbit Integrations 
3.3.1. Methodology & Overview 

To get a better understanding of the past orbital his- 
tory of Leo I we need to use more detailed models for 
the MWs gravitational potential $mw- Following Besla 
et al. (2007), we describe this potential as a static, ax- 
isymmetric, three-component model consisting of dark 
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Table 2 

Model Parameters for the Milky Way 



M vir a 




R v i T C 


M disk d 


(M e ) 


Cvir 


(kpc) 


(Af©) 


1.0 x 10 12 


9.86 


261 


6.5 x 10 10 


1.5 x 10 12 


9.56 


299 


5.5 x 10 10 


2.0 x 10 12 


9.36 


329 


5.0 x 10 10 



a Mass contained within the virial ra- 
dius. 

b The virial concentration parameter 
(Klypin et al. 2011). 

c The virial radius. See text for defini- 
tion. 

d Mass of the disks. 

matter (DM) halo, disk (Miyamoto & Nagai 1975) and 
stellar bulge (Hcrnquist 1990): 

*MW = $DM halo + 'I'disk + ^bulgo- (10) 

The DM halo is initially modeled as an NFW halo 
(Navarro et al. 1997) with a virial concentration parame- 
ter (c v ; r ) defined as in Klypin et al. (2011) from the Bol- 
shoi Simulation (see also, van der Marel et al. 2012b). 
We apply the adiabatic contraction of the NFW halo in 
response to the slow growth of an exponential disk using 
CONTRA code (Gnedin et al. 2004). The density profile 
of the MW is then truncated at the virial radius 8 . The 
bulge mass of the MW is kept fixed at 1O 1O M with a 
Hernquist scale radius of 0.7 kpc. The exponential disk 
scale length is also kept fixed at 3.5 kpc. 

We adopt three different mass models for the MW with 
total virial masses of 1.0 x 10 12 , 1.5 x 10 12 , and 2.0 x 
10 12 Mq. In all cases, the bulge is modeled with a scale 
radius of 0.7 kpc and a total mass of 10 10 Mq. The disk 
scale radius is also kept fixed at 3.5 kpc, but the disk mass 
is allowed to vary to reproduce the circular velocity at 
the Solar circle. We adopt a MW circular velocity of 239 
km s" 1 at the solar radius of 8.29 kpc (McMillan 2011). 
Model parameters are listed in Table 2. 

The escape velocities at the distance r = 260.6 kpc of 
Leo I are 182, 222, and 256 kms^ 1 for the models with 
masses of 1.0 x 10 12 , 1.5 x 10 12 , and 2.0 x 1O 12 M , re- 
spectively. 9 For the latter two models, this exceeds 
our best estimate v — 196.0 ± 19.4 km s -1 for the total 
Galactocentric velocity of Leo I. So Leo I is most likely 
bound to the MW if M M w,vir > 1-5 x 10 12 M Q , as al- 
ready suggested in Section 3.2. For the lowest-mass MW 
model studied here, with MM\v,vir = 1.0 x 10 12 M Q , Leo I 
is on an unbound hyperbolic orbit. This has repercus- 
sions for the viability of such low mass MW models in 
a cosmological context, since satellites are rarely found 
on hyperbolic orbits. We explore this in more detail in 
Paper II. 

Using the current Galactocentric position and velocity 

8 The virial radius is denned as the radius such that p v ; r = 
A v i r Q m p cr i t , where the average overdensity A v i r = 360 and the 
mean matter density parameter S7 m = 0.27 (see equation Al in 
van der Marel et al. 2012a). 

9 In a halo that is not truncated at the virial radius, the escape 
velocities are larger by 45-55 kms -1 . At fixed Mmw,vii, this in- 
creases the probability that Leo I is bound; or at fixed probability, 
this means that Leo I is bound even at lower Mmw vir- While 
the escape velocity of an NFW halo is finite, its mass is not. So 
truncation at some large radius is always physically motivated. 



vectors of Leo I as initial conditions, we can solve the 
differential equations of motion numerically to follow the 
velocity and position of Leo I backward in time. If we 
consider only the gravitational influence of the MW, the 
equation of motion has the form: 

* = ^*mw(M), (11) 

In Section 3.3.2 this equation of motion is solved using 
well-established numerical methods in order to constrain 
Leo Fs interaction history with the MW (i.e., pericen- 
tric distance and epoch of accretion) over a Hubble time. 
Leo I is modeled as a Plummer potential with a soften- 
ing length of 0.5 kpc, and total mass of 1.3 x 1O 8 M (see 
Table 3). With these parameters, the dynamical mass of 
Leo I within 0.93 kpc is 8.9 x 10 7 M Q , as expected from 
Table 2 of Walker et al. (2009) (i.e. the inferred mass 
within the outermost data point of the empirical veloc- 
ity dispersion profile, referred to as ri ast in Walker et al. 
2009). We note that dynamical friction is expected to 
be negligible for such a low mass satellite (even when its 
extended dark mass outside the outermost data point is 
included) and it is thus not computed in the equation of 
motion. 

Other members in the Local Group with significant 
masses may exert dynamical influence on the orbital his- 
tory of Leo I. Relevant to our analysis are the LMC and 
M31. It has been theorized that a number of satellite 
galaxies of the MW lie in a similar orbital plane as the 
Magellanic Clouds, referred to as the "Magellanic Plane" 
of galaxies (Kunkel & Demers 1976; Lynden-Bell 1976). 
The ubiquity of this statement is directly relevant to cos- 
mological studies of how satellite galaxies are accreted by 
MW mass halos (e.g., D'Onghia & Lake 2008; Metz et al. 
2009; Sales et al. 2011; Wang et al. 2012). The Sagit- 
tarius dSph is known to be orbiting in a plane that is 
perpendicular to that of the Magellanic Clouds (which 
orbit approximately in the Y — Z Galactocentric plane). 
As such, it is clear that at least one of the MW satellites 
is an outlier. Without accurate proper motions, it has 
been unclear to what extent Leo Fs orbit lies in a different 
plane or has a different rotational sense. Furthermore, 
the Keplerian orbit analysis of Section 3.2 demonstrates 
that Leo I is likely on a fairly eccentric orbit. On such a 
high eccentricity orbit, it is possible that M31 may exert 
an important gravitational influence at early times. 

In Section 3.3.4 we compute the orbital histories of the 
LMC and M31 in addition to that of Leo I to define its 
orbital plane relative to that of the LMC and also to 
assess the dynamical significance of M31 to its orbital 
history and origin. We modify the equations of motion 
(Eq. 11) to include the gravitational influence of both 
M31 and the LMC while computing the orbital history 
of Leo I over the past 8 Gyr. We simultaneously compute 
the equations of motion for each of the MW, M31 and 
LMC, accounting for the gravitational acceleration from 
the other bodies. The orbits are computed in the Galac- 
tocentric frame, using the current velocities and positions 
of the LMC (Kallivayalil et al. 2012) and M31 (van der 
Marel et al. 2012a), as listed in Table 3. The masses 
of the MW and M31 are assumed to be static over time, 
which limits the accuracy of this analysis past 6 Gyr (i.e., 
when the mass of these galaxies is expected to be about 
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Table 3 

Properties of Local Group Galaxies 





M,„ t ,i 

i,x total 




r last b 


M(< ri aat ) c 

\ ^ last ) 


X d 


yd 


Z d 


V x e 


vy° 


V z e 


Ref. f 




Ha 1 avv 
v. i <x,i (i - k v 


(Mr-,) 


(kpc) (kpc) 


(Mr-,) 






(kne) 


(km s _1 1 (k 


m s _1 1 Ck- 

Lll tD ) \^-- 




Dist. 


RV 


PM 


Leo I 


1.30 x 10 8 


0.5 


0.93 


8.9 x 10 7 


-125.0 


-120.8 


194.1 


-167.7 


-37.0 


94.4 


1,2,3 


4 


This study 


Carina 


5.70 x 10 7 


0.5 


0.87 


3.7 x 10 7 


-24.8 


-94.8 


-39.3 


-72.9 


6.9 


38.0 


2,5 


6 


7 


Draco 


3.90 x 10 8 


0.5 


0.92 


2.6 x 10 8 


-3.5 


76.3 


53.0 


17.1 


56.0 - 


-227.9 


8 


9 


10 


Fornax 


1.45 x 10 8 


0.5 


1.70 


1.3 x 10 8 


-40.0 


-49.2 


-129.4 


-24.5 


-140.7 


106.2 


5,11,12,13,14,15 


16 


17 


Leo II 


6.40 x 10 7 


0.5 


0.42 


1.7 x 10 7 


-77.3 


-58.3 


215.3 


102.2 


237.0 


118.4 


18,19 


20 


21 


Sculptor 


1.35 x 10 8 


0.5 


1.10 


1.0 x 10 8 


-5.3 


-9.6 


-84.1 


-19.4 


224.6 - 


101.6 


13,22 


23 


24 


Sextans 


2.90 x 10 7 


0.5 


1.00 


2.0 x 10 7 


-40.0 


-63.6 


64.6 


-181.1 


113.6 


113.6 


25 


26 


27 


Ursa Minor 


7.70 x 10 7 


0.5 


0.74 


4.4 x 10 7 


-22.2 


52.1 


53.6 


-107.5 


-15.2 


116.1 


8 


9 


28 


LMC 


5.00 x 10 10 


11.0 


9.00 


1.3 x 10 10 


-1.1 


-41.0 


-27.8 


-57.4 - 


-225.6 


220.7 


29 


30 


31 


M31 


See text 








-378.9 


612.7 


-283.1 


66.1 


-76.3 


45.1 


32 


33 


34 



a Softening parameter for the Plummer profile. 

b Radius of the outermost data point of the empirical velocity dispersion profile for each satellite, as defined in Walker et al. (2009). 

c Mass inferred within ri ast . Note that for the LMC, this is the mass within the radius of the last data point in the carbon star analysis of 

van der Marel et al. (2002). 

d Galactocentric positions with the origin at the Galactic center, the Z-axis pointing toward the Galactic north pole, the X-axis pointing in 
the direction from the Sun to the Galactic center, and the Y-axis pointing in the direction of the Sun's Galactic rotation. 
c Galactocentric velocities with vectors pointing toward X, Y, and Z as defined above. 

1 Data for distance (Dist.), heliocentric radial velocities (RV), and proper motions (PM) were taken from the following references. (1) Mendez 
et al. 2002; (2) Bellazzini et al. 2004; (3) Held et al. 2010; (4) Mateo et al. 2008; (5) Pietrzynski et al. 2009; (6) Mateo et al. 1993; (7) Piatek 
et al. 2003; (8) Bellazzini et al. 2002; (9) Armandroff et al. 1995; (10) Scholz & Irwin 1994; (11) Bersier 2000; (12) Saviane et al. 2000; (13) 
Rizzi et al. 2007a; (14) Rizzi et al. 2007b; (15) Gullieuszik et al. 2007; (16) Walker et al. 2006; (17) Piatek et al. 2007; (18) Bellazzini et al. 
2005; (19) Gullieuszik et al. 2008; (20) Koch et al. 2007; (21) Lepine et al. 2011; (22) Pietrzynski et al. 2008; (23) Qucloz et al. 1995; (24) 
Piatek et al. 2006; (25) Lee et al. 2003; (26) Hargreaves et al. 1994; (27) Walker et al. 2008; (28) Piatek et al. 2005; (29) Freedman et al. 
2001; (30) van der Marel et al. 2002; (31) Kallivayalil et al. 2012; (32) Freedman & Madore 1990; (33) Courteau & van den Bergh 1999; (34) 
Sohn et al. 2012. 



half of their current value, e.g. Fakhouri et al. 2010). 

M31 is modeled using an NFW halo, where its virial 
mass is determined by preserving the total mass of the 
Local Group, given the mass of the MW in each model. 
The density profile of M31 is also modeled to be trun- 
cated at the virial radius. Using the recent proper mo- 
tions of M31 by Sohn et al. (2012) and other mass ar- 
guments in the literature, van der Marel et al. (2012a) 
estimate the Local Group mass to be Mlg = (3.17 ± 
0.57) x 1O 12 M . In this analysis we thus require that the 
combined Mmw + Mm3i = 3 x 10 12 M©. Given the large 
distances involved, the contribution of M31's disk/bulge 
component to the gravitational influence on Leo I is ir- 
relevant. 

The LMC is modeled as a Plummer sphere with a soft- 
ening parameter of 11 kpc and total mass of 5 x 1O 1O M . 
With these parameters, the mass contained within 9 kpc 
is <~ 1.3 x 1O 1O M0, as observed (van der Marel et al. 
2002). To model the orbital evolution of a massive satel- 
lite such as the LMC accurately, dynamical friction ef- 
fects owing to its motion through the dark matter halo 
of the MW must be accounted for. Dynamical friction is 
included using the Chandrasckhar formula with an ap- 
proximation to the Coulomb logarithm as in Besla et al. 
(2007). Meanwhile, dynamical friction is irrelevant to 
the motion of the MW and M31 because their halos do 
not overlap over the past 8 Gyr. 

The star formation history of Leo I has been the topic 
of many studies, especially with HST (Caputo et al. 
1999; Gallart 1999; Dolphin 2002; Smecker-Hane et al. 
2010). From the most recent HST ACS/WFC obser- 
vations (Smecker-Hane et al. 2010), Leo I is known to 
have formed stars continuously since > 12 Gyr ago, with 
two pronounced star formation activities at ~ 4.5 and 
<~ 2 Gyrs ago. After this last activity, star formation 



abruptly dropped until a complete cessation at ~ 0.5 
Gyr ago. Some of the inferred increases and decreases 
in star formation activity may be related to features in 
Leo Fs orbit about the MW, including the time of accre- 
tion and the pericenter time. We determine these times 
in Section 3.3.2. 

The origin of enhanced star formation activities in 
Leo Fs past may also be related to interaction with other 
satellites. Furthermore, three-body encounters may also 
alter the orbital trajectory of Leo I as proposed by Mateo 
et al. (2008), potentially explaining its high speed today 
(e.g., Sales et al. 2007). To test for possible interac- 
tions with other MW satellites, we extend in Section 3.3.5 
the analysis of Section 3.3.2 such that the equations of 
motion now account for gravitational interaction terms 
with the following satellites for which proper motions are 
available from other studies (see Table 3): Carina, Draco, 
Fornax, Leo II, Sculptor, Sextans, Ursa Minor, and the 
LMC. The Sagittarius dSph is not included in this analy- 
sis because its orbit is too close to the MW disk plane to 
dynamically influence that of Leo I. The SMC is also not 
included because its orbit is likely closely matched to that 
of its binary companion, the LAIC. Since the LMC is the 
more massive of the pair, it is likely to be the more signif- 
icant perturber. Each satellite is modeled as a Plummer 
potential with total mass and softening parameters as 
listed in Table 3; model parameters are chosen to match 
the observed masses within ri as t , as defined in Walker et 
al. (2009). 

3.3.2. Leo I Orbital Properties 

Here we present the plausible orbital histories of Leo I 
following the methodology outlined above. To propagate 
the observed errors, we randomly sampled 10,000 combi- 
nations of the west and north components of the observed 
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Figure 6. Ten thousand points randomly sampled from the 4<r proper-motion error space of Leo I. The dashed ellipses indicate the 
standard deviation of the enclosed points from the mean (filled triangle). For each point, the orbital history of Leo I was computed by 
integrating the equations of motion [equation (11)] backward in time for the lowest mass MW model (M v „ = 1.0 X 10 12 M Q ). In panel (a), 
solutions with more than one previous pericentric approach are highlighted. All cases have at least one recent pericentric approach with 
the MW. In panel (b), points are color coded based on the distance of their most recent pericentric approach, as indicated in the legend. 
The minimum (maximum) pericentric approach is indicated by the black open square (pentagon). In panel (c), the color coding indicates 
the time of the pericentric approach. Finally, panel (d) indicates the infall time, i.e. when Leo I last entered the virial radius of the MW. 



Leo I proper motion within the error space provided in 
Section 2.3. For each of these combinations, the 3D ve- 
locity was derived and the orbit of Leo I was followed 
backward in time for a Hubble time in each of our three 
MW models. The analysis presented in this section only 
considers the influence of the MW's gravitational field on 
the orbit of Leo I, as described in equation (11). 

We seek to define Leo I's interaction history with the 
MW by constraining the time and distance of Leo I's re- 
cent pericentric approach to our MW (r pcr i, t peT i), the 
number of such encounters it might have had in the past 
(-Wperi) and the time of infall to our system (Onfall)- We 
define infall time as the time at which Leo I last crossed 
the virial radius of the MW. The allowed range in peri- 
center distance will inform us about the maximal tidal 



influence the MW may have exerted over Leo I, which 
may have caused its transformation into a dSph or trig- 
gered enhanced star formation activities. The pericenter 
distance also informs us about the importance of ram 
pressure stripping by the MW's gaseous halo in removing 
gas from this system (e.g., Grcevich & Putman 2009); the 
deeper that Leo I travels into the MW's halo, the higher 
the background gas density and the more likely gas gets 
stripped. The infall time is similarly relevant; the longer 
ago that Leo I was accreted, the longer the time scale for 
ram pressure stripping to operate. Ultimately, the num- 
ber of previous pericentric passages will tell us whether 
Leo I is in fact a recent interloper in our system. These 
properties can be further used to constrain the halo mass 
of our own MW statistically by identifying analogs of the 
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Leo I satellite (in terms of mass and orbital properties) 
about MW-type hosts in large-scale cosmological simu- 
lations (see Paper II). 

Figures 6, 7 and 8 show the 4a proper motion error 
space that is sampled to determine the orbital properties 
of Leo I in the three different MW mass models. Points 
are color coded to reflect the range of allowed values for 
the quantity of interest (7Vp er i, r pe ri, ^peri and ijnfaii 

). The 

mean values for these quantities are listed in Table 4. 
In addition, we also list the mean velocity at pericenter 
relative to the escape speed at r per i, (v pei i/v esc ), and the 
velocity at infall relative to the circular velocity of the 
MW at the virial radius (winfall/^Rvir)' In the cases where 
there is a second pericentric approach some time in the 
past, we list the time this occurs (i per i,2) and the mean 
apocenter distance (f a po)- We mark the minimum and 
maximum rp er i allowed within the 4<r proper motion error 
space and the times at which they occur in the second 
and third panels in each of Figures 6, 7 and 8 by the 
black open square (min) and pentagon (max). 



In all cases, Leo I has recently had a pericentric passage 
with respect to the MW, and so panels (a) in Figures 6, 7 
and 8 only note if a second pericentric approach occurs. 
As the MW mass increases, the likelihood of a second 
close passage also increases; however there are no cases 
where a third pericentric approach occurs. In the lowest 
MW mass model, there are solutions for a second pericen- 
tric passage only outside the 2a error ellipse. However, 
there are only two such solutions out of our 10,000 re- 
alizations, one with a second pericenter time of 13 Gyr 
and the other with 2 Gyr (which is likely a slingshot 
scenario where Leo I got too close to the MW center). 
This reflects the fact that for the low-mass MW model, 
Leo I is generally on a hyperbolic or near-hyperbolic or- 
bit. Solutions for a second pericenter are readily ob- 
tainable within the la error ellipse for the higher MW 
mass models. However, the Monte-Carlo statistics (see 
Table 4) still favor orbits with only one previous pericen- 
ter. Moreover, in those cases with a second pericenter, 
the time since that pericenter is > 10 Gyr. Also, the 
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Figure 8. Similar to Figure 6, but now for the high mass MW model (Af v ; r = 2.0 X 10 12 ). 



Table 4 

Mean Orbital Properties of Leo I 



Mmw 
(xlO 12 M ) 


(-^peri) 


{''peri) 

(kpc) 


{^pcri ) 

(Gyr) 


(^peri/^esc) 


(ipo> 

(kpc) 


(*peri,2) 

(Gyr) 


(Onfall) 

(Gyr) 


(I'mfall/fRvir} 


1.0 
1.5 
2.0 


1.0 ± 0.2 
1.0 ±0.2 
1.5 ±0.5 


100 ± 37 
90 ±35 
82 ±33 


1.1 ± 0.1 
1.1 ±0.1 
1.0±0.1 


1.04 ±0.07 
0.96 ±0.04 
0.91 ±0.02 


1177 ± 640 
558 ± 61 
530 ± 55 


7.7 ± 5.6 
12.6 ± 1.3 
11.0 ±1.7 


2.2 ±0.2 

2.3 ±0.1 
2.5 ±0.2 


1.6 ±0.2 
1.3 ±0.2 
1.0 ±0.1 



Note. — Orbital properties as defined in the text. Standard deviations for the mean values arc marked. Apocenter 
distances are computed only for the cases that have a second pericentric approach. 



apocenter distances are well beyond the virial radius of 
the MW for all models (d apo > 500 kpc). Recall that 
the MW mass was assumed to be static in time; with 
an accurate treatment of the mass evolution of the MW 
it is doubtful that any of these second pericentric pas- 
sages would still occur. It is thus most likely that Leo I 
has passed its first infall into the MW. Moreover, if the 
previous pericentric passage occurred > 10 Gyr ago, this 



implies that Leo I would be at large distances from the 
MW exactly a Hubble time ago, which is an implausible 
scenario in the view of the timing argument. 

The average time of Leo I's most recent pericentric 
passage is ~1 Gyr ago for all MW models. This is 
similar to what was found from the Keplerian calcula- 
tions. The average pericentric approach ranges from 80- 
100 kpc, and the average velocity at pericenter ranges 
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from 300-370 km s _1 (Table 4). This distance is some- 
what larger, and the velocity somewhat smaller, than 
found from the Keplerian calculations. This is easy to 
understand from the fact that a Keplerian model is too 
concentrated, and therefore overestimates the accelera- 
tion as Leo I approaches the MW. The average ratio of 
pericenter to apocenter distance in the orbit calculations 
ranges from 0.08-0.15. This is larger than in the Ke- 
plerian calculations, in part because those calculations 
ignore the mass outside of the Leo I distance, and there- 
fore overestimate the apocenter distance. 

The tidal radius r t = r po ri[2M M w(r P eri) /-Mlco i]~ 1/3 
of Leo I, given pericentric distances of 80-100 kpc, is 
3-4 kpc for the three MW models. The present optical 
radius of Leo I is ~ 1 kpc (Sohn et al. 2007; Walker et 
al. 2009). This means that on average, the tidal field 
of the MW does not appear to be sufficient to tidally 
truncate Leo I to this radius, implying that Leo I likely 
has an extended dark matter halo. There is a clear trend 
in r pcr i with the west component of the proper motion 
as evidenced in panel (b) of Figures 6, 7 and 8; r por i 
increases with increasing fiw- The minimum pericentric 
approach determined in the 10,000 Monte-Carlo orbits is 
~1— 2 kpc. However, a pericentric approach of < 5 kpc 
is unlikely, because the tidal radius of Leo I in all MW 
mass models is less that 0.7 kpc. Leo I is not sufficiently 
disturbed to have approached this close to the MW. On 
the other hand, pericentric approaches of 20 kpc yield 
tidal radii of ~1 kpc. Cases where 15 < r pcr i < 25 kpc 
are found within the 1.5-2tr error ellipse in all MW mass 
models. So such close approaches are not ruled out by our 
data. However, they are not so likely given our data, with 
only 2-3% of the Monte-Carlo calculated orbits yielding 
r pori < 25 kpc. 

The time of Leo I's last pericentric passage, ~ 1 Gyr 
ago, corresponds roughly to the time when star forma- 
tion stopped in Leo I (Caputo et al. 1999; Gallart 1999; 
Dolphin 2002; Smecker-Hane et al. 2010). The pericen- 
tric approach is the point in time where Leo I would 
experience maximal ram pressure stripping, which could 
lead to quenching. All satellites of the MW within 300 
kpc (apart from the Magellanic Clouds) are devoid of gas 
(Grcevich & Putman 2009), consistent with this picture. 
As the MW mass increases, i per i decreases; for the most 
massive MW model the maximal i per i is 1.3 Gyr and the 
minimal value is 0.6 Gyr. These values are remarkably 
consistent regardless of MW mass, when searching the 
full 4cr proper motion error space [panels (c)]. Hence, 
it is likely that star formation stopped in Leo I owing 
to ram pressure effects at pericentric approach, implying 
a quenching time scale of ii n f a ii — t pcll rj 1.3 Gyr (see 
Table 4). 

Of course, in general, star formation can cease in galax- 
ies for many other reasons, e.g., exhaustion or blowout 
of the gas supply. However, if this were the cause of 
star formation ceasing in Leo I, then there would be no 
natural explanation for why this would coincide with a 
pericenter passage. On the other hand, this coincidence 
could certainly happen by chance, especially since the 
uncertainties in both the SFH and the orbital analysis 
are significant. 

The average infall time ranges from 2.2-2.5 Gyr with 
little variation, regardless of MW mass. The infall time 
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Figure 9. Comparison between the average Leo I proper motion 
measured in this study (red data point with error bars) and the 
proper motions predicted by the models of Sohn et al. (2007) (green 
closed square) and Mateo et al. (2008) (blue open diamonds). The 
black cross indicates the proper motion for a radial orbit, i.e., with 
Vtan = 0, as in Figure 5. 

is similar to the time scale of the most recent enhanced 
star formation observed in Leo I (Caputo et al. 1999; 
Gallart 1999; Dolphin 2002; Smecker-Hane et al. 2010), 
suggesting that this enhanced star formation activity was 
triggered by either ram pressure compression as Leo I 
entered a higher gas density environment relative to the 
Local Group, or as it began to feel gravitational torques 
exerted by the MW. Note that imfaii refers to the most 
recent time at which Leo I entered the virial radius; there 
are cases where Leo I has made an earlier pericentric 
approach about the MW. But, as discussed earlier, such 
orbits may not be physical or plausible. There are a few 
cases where Leo I remains within the virial radius of the 
MW for approximately a Hubble Time [indicated by red 
squares in panels (d)]. However, such a scenario has low 
likelihood, since these cases are all 4cr outliers that only 
occur in the high mass MW model. 

The ratio between the infall velocity and the circular 
velocity at the virial radius (i>R V ir) ranges from 1.0-1.6 
(Table 4). For the low mass MW model, the average 
Leo I infall velocities are higher than expected based on 
cosmological simulations of structure formation, where 
subhalos are typically accreted with characteristic or- 
bital velocities of ~ l.li>R V i r at the virial radius (lcr scat- 
ter of 25%, Wetzel 2011). This further disfavors masses 
AWvir < 10 12 Af Q . 

Our conclusion that Leo I has most likely passed its 
first infall into the MW, and our value for ii n f a u are con- 
sistent with the results of Rocha et al. (2012) for Leo I. 
They find that there is a tight correlation between the 
present day orbital energies of the MW satellites and 
their infall times as inferred from cosmological simula- 
tions. We explore the implications of this for our under- 
standing of Leo I further in Paper II. 

3.3.3. Comparison to Previous Orbit Estimates 
Sohn et al. (2007) and Mateo et al. (2008) provided es- 
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Figure 10. The mean orbital history of Leo I for the low mass MW model (M v j r = 1.0 X 1O 12 M0). The mass of M31 is chosen to preserve 
the total Local Group mass of 3 X 10 12 Mq. Panel (a) shows the separation between Leo I and the LMC, MW and M31 as a function of 
time. In panels (b), (c), and (d), the orbital plane is presented in the Galactocentric X — Y , Y — Z , and X — Z planes, respectively. The 
virial radius (R v ir) of the MW is indicated by the dashed horizontal line in panel (a), and the dot-dashed ellipses in the other panels. In 
panels (b)-(d), the MW location is indicated by the thick black line (circle in the case of the Y — Z plane) centered at (0,0). The current 
location of each galaxy is indicated by the solid squares. Solid black lines indicate results when Leo I is assumed to be moving with the 
mean proper motion determined in this study. For comparison, we also plot the orbit of Leo I accounting only for the influence of the MW 
in dashed green lines. The thin dashed-dotted black lines indicate solutions where the pericenter distance is minimized and maximized [3<r 
outliers; see panel (b) in Figure 6]. 



timates of the orbital history of Leo I based on indirect 
arguments, rather than proper motion measurements. 
They aimed in particular to explain the photometric and 
kinematical data of giants stars in Leo I. The proper mo- 
tions corresponding to the proposed orbits are compared 
to our new HST measurements in Figure 9. 10 The or- 
bits of Sohn et al. (2007) had a pericentric approach of 
only 10 kpc, so the predicted proper motion is very 

10 The proper motion predictions reported in Sohn et al. (2007) 
are erroneous. We re-derived the predicted proper motions based 
on the orbital positions and velocities of their model 117 at t = 
and the result is (/iw,/in) = (0.032,-0.160) masyr - 1 . This 
assumed the same position and velocity of the Sun as the present 
study. 



close to the Vtan = point. Mateo et al. (2008) provided 
proper motion predictions for a range of assumed Leo I 
masses (their Table 8). While their predicted perigalac- 
tic distances reach out to 30 kpc, their proper motions lie 
on the opposite side of the Vt an = point compared to 
our measurements. So our measurements do not confirm 
the predictions of these studies. More specifically, the 
previous studies argued for highly eccentric orbits with 
smaller perigalactic distances than what we find here. 

Sohn et al. (2007) focused primarily on trying to repro- 
duce the observed photometric and kinematic features 
of Leo I by adopting a tidal disruption scenario. The 
orientation of their model orbits was determined by as- 
suming that the position angle of the Leo I ellipticity 
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Figure 11. Similar to Figure 10, but now for the intermediate MW mass model (M v ; r = 1.5 X 1O 12 M0). 



and the orientation of the break population are caused 
by tidal effects and tidal stripping, respectively. They 
showed that the observed features can be plausibly pro- 
duced by the tidal effects of the MW. However, the tidal 
effects may have been overestimated, given that the or- 
bital properties derived here imply that Leo I is on a less 
eccentric orbit than assumed by Sohn et al. (2007). This 
discrepancy does not imply though that the tidal scenar- 
ios used by Sohn et al. (2007) and Mateo et al. (2008) 
are necessarily wrong. It may just be that some of the 
specific assumptions in their models were oversimplified. 
For example, if Sohn et al. (2007) had not modeled the 
progenitor satellite as a spherical Plummer profile, the 
best-fit orbits may well have been more consistent with 
the observed proper motion. New iV-body models based 
on the observed proper motion should be able to fur- 
ther improve our understanding of the tidal disruption 
features observed in Leo I. However, such models are be- 
yond the scope of the present paper. 



3.3.4. Leo I Orbital Plane 

Here we compute the orbital history of Leo I, relative 
to the other major players in the Local Group, namely 
the MW, the LMC, and M31. We aim to define Leo I's 
orbital plane and compare it to that of the LMC, and to 
determine whether Leo I was ever close enough to M31 
for it to have exerted any dynamical influence in Leo I's 
history. 

In Figures 10, 11, and 12 we plot the orbit of Leo I 
using the mean proper motions determined in this study, 
including the influence of the MW, LMC, and M31 (solid 
black line). For comparison, we also plot the orbit 
of Leo I accounting only for the influence of the MW 
(dashed green line). Orbits corresponding to proper mo- 
tions that are ±3er from the mean (identified as having 
min/max pericenter distances to the MW from panels (b) 
in Figures 6, 7, and 8) are indicated by the thin dash- 
dotted black lines (Rmin and Rmax) in all panels. 

Panels (a) illustrate the separation of Leo I from the 
other galaxies and Panels (b), (c), and (d) respectively 
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Figure 12. Same to Figure 10, but now for the high MW mass model (M vir = 2.0 X 1O 12 M ). 



show the orbits of the galaxies in the X-Y, Y-Z, and 
X-Z Galactocentric planes. As the MW mass increases, 
Leo Fs past orbit becomes less eccentric and increasingly 
directed towards M31. It is clear that Leo I does not get 
closer than 400 kpc from M31 in any model over the past 
8 Gyr and neither the presence of M31 nor the LMC have 
an impact on the infall time or pericenter properties 11 . 
However, M31 may play a role at early times in the higher 
MW mass models; the orbits are more energetic, reach- 
ing larger distances than if the gravity of the MW were 
considered alone. We note that we have not accounted 
for the mass evolution of the MW or M31, which will 
diminish the role that M31 plays at early times. But, 
this analysis does suggest that accounting for the local 
overdensity (i.e. that there are two roughly equal mass 
galaxies in our Local Group) may be a relevant parameter 
in understanding the origin of the angular momentum of 

11 This is unsurprising given that these properties are computed 
using a backwards integration scheme. 



high speed satellites. In this study, such considerations 
are only relevant for the higher mass MW models; in 
the lowest mass MW model, Leo I is too far from M31 
at all times for torques to be relevant. Our conclusions 
regarding the hyperbolic nature of Leo I's orbit in low- 
mass MW models is thus robust to the influence of M31 
and can be compared statistically to satellite orbits found 
in cosmological simulations of isolated MW analogs (see 
Paper II). 

The minimum pericentric passage orbits appear to be 
slingshot orbits, approaching near the MW center, gain- 
ing energy and traveling to larger distances. However, as 
discussed in Section 3.3.2, such orbits are likely unphys- 
ical because Leo I is not sufficiently disturbed to have 
traveled this close to the MW. Since the orbit of Leo I is 
still bound to the MW in the vast majority of models, we 
do not need a small pericenter, slingshot orbit to explain 
the orbital properties of Leo I. 

Panels (b), (c), and (d) respectively show the orbits of 
the galaxies in the X-Y, Y-Z, and X-Z Galactocentric 
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Figure 13. The minimum separation to Leo I that is allowed within the proper motion error space of each satellite is plotted as a function 
of time for each MW mass model. Cases where the separation is sufficiently small so as to tidally truncate Leo I to within at least 2 kpc 
are marked in blue. Given the low mass ratios, such separations will also allow Leo I to influence that satellite in turn, except for the 
LMC. Black dashed lines indicate known epochs of enhanced star formation activity in Leo I, as discussed in the text. The red dashed line 
indicates its infall time within the virial radius of the MW. The size of each point indicates the distance of the corresponding Leo I proper 
motion from the mean observed proper motion, as described in the legend; smaller points are more likely. Although events were found 
within ltr for some satellites, they typically represent only < 4% of the Monte-Carlo drawings. This is shown explicitly in Figure 14 for 
the satellites that have the closest encounters, Ursa Minor and Leo II, to illustrate the low probability of such encounters. We note that 
Sculptor never gets closer to Leo I than its current separation and is thus omitted from the plot. 



planes. The orbital angular momentum of Leo I is not 
coincident with that of the LMC. This is most clearly 
illustrated when looking at the orbital history in the X-Z 
Galactocentric plane, especially in the higher-mass MW 
models. 12 The LMC is moving in a clockwise direction 
in this plane whereas Leo I is moving counterclockwise. 

3.3.5. Interactions with other Satellites 

We have so far established that Leo I is likely on its first 
orbit around the MW and that its most recent pericen- 
tric approach was likely at a separation of 80-100 kpc, 
too large for the MW to have exerted significant tidal 
torques. Yet, Leo I shows signs of a past interaction. 
Sohn et al. (2007) found an excess of red giant stars along 
the major axis of Leo I's main body relative to a symmet- 
ric King profile. In addition to this spatial configuration, 
Leo I red giant stars have an asymmetric radial velocity 
distribution at large radii (cf., see also Mateo et al. 2008). 
If the MW is not the culprit for the distorted structure 
and kinematics of Leo I, then what is? 

Here we consider the orbital history of not only Leo I, 
but also of the other satellite galaxies of the MW simul- 
taneously. We randomly sampled 10,000 combinations of 
the observed west and north components of the proper 
motion within the 4<r error space for Leo I and for each 
satellite. The proper motions, distances, and radial ve- 
locities were taken from various sources as listed in Ta- 
ble 3. As for Leo I, we took the distance for each satellite 
to be the error- weighted average of TRGB measurements 
in the last decade. The proper motions and radial veloc- 
ities were adopted from the most recent measurements 
available in the literature. We computed the orbits of 
all satellites backwards in time for 10 Gyr using each of 
our three MW models, and using the mass model for each 
satellite outlined in Table 3. Our goal is to determine the 

12 In the lower-mass MW models (Figure lOd) the LMC and 
Leo I are also clearly unassociated, because they are more than 
300 kpc from each other for most of their evolution. However, the 
LMC orbit has a less clear sense of rotation in the X-Z plane, owing 
to the gravitational influence of M31 which causes a kink/twist in 
the orbit. The mass of M31 is the highest in the lowest MW mass 
model, thereby making its gravitational perturbation the strongest. 



closest separation that Leo I may have reached to any of 
these other satellites within the error space and to assess 
whether such separations are sufficient to exert torques 
on Leo I and induce star formation or to significantly 
modify Leo I's orbit. 

Given the small masses of these satellites, their dynam- 
ical influence is minimal unless the separation between 
them is small. Overall, satellite separations as low as 3- 
4 kpc are required before one can significantly influence 
the other, i.e. such that the tidal radius of the satellite 
is less than 2 kpc. Because the LMC is much more mas- 
sive than the other satellites, a larger separation of at 
least 20 kpc will allow it to distort Leo I to within 2 kpc. 
Of course, Leo I is much too small to strongly affect the 
LMC 13 . 

Figure 13 indicates the minimum separation between 
each satellite and Leo I as a function of time for each 
respective MW mass model. Satellites that reach sepa- 
rations within their 4cr error space small enough to in- 
fluence Leo I within a radius of 2 kpc (and vice versa) 
are highlighted in blue. The size of the point reflects the 
probability of that encounter, based on the sigma devia- 
tion of the required Leo I proper motion from the mean. 
Vertical dashed lines indicate relevant events in the his- 
tory of Leo I, such as epochs of star formation and the 
epoch of infall into the MW. 

There are no obvious satellite encounters found that 
can explain the enhanced star formation activities at ~ 2 
and ~ 4.5 Gyr ago. It appears much more likely that 
the enhancement of star formation that took place ~ 2 
Gyr ago is related to the epoch of accretion by the MW. 
While there are some parameter combinations that allow 
for close encounters between Leo I and Fornax, Carina 
and the LMC, these are likely random events as there are 
only ever 2, or 3 points within the 10,000 orbit search 
that yield such encounters. Also, Sculptor never gets 
close to Leo I than it is today, and is thus omitted from 
the plot. Among the MW satellites we consider, there 
is a higher probability - though still relatively small - 

13 Indeed the SMC's tidal field has had limited influence over 
the star formation history of the LMC (Besla et al. 2012). 
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M„ (mas/yr) M w (mas/yr) 

Figure 14. Orbits, out of 10,000 randomly selected points from the 4cr error space of the Leo I proper motions, that produce a passage 
of Ursa Minor (left panel) or Leo II (right panel) within 10 kpc of Leo I. Orbits are color coded based on the minimum separation, as 
calculated for the intermediate MW mass model (M v i r = 1.5 X 10 12 Af(j). Separations less than 3 kpc are readily found within the lcr error 
ellipse, indicating that such past encounters are possible. For Ursa Minor, there are 62 blue points and 19 red points that fall within the 
1<t region (the innermost circle), and for Leo II, there are 171 blues points and 57 red points within the la region. 



of an encounter with Ursa Minor and Leo II for all MW 
models. Specifically, 1.3% and 4.5% of the Monte-Carlo 
calculated orbits yield passages within 10 kpc for Ursa 
Minor and Leo II, respectively. These encounters are ex- 
pected to have occurred within the past 1 Gyr, well after 
star formation has ceased in all of these galaxies (Carrera 
et al. 2002; Koch et al. 2007; Kirby et al. 2011). The en- 
counter is thus unlikely to have signatures in the star for- 
mation histories of these galaxies, as they should already 
have been devoid of gas at that point. However, there 
may be kinematical signatures instead. In particular the 
distortions in Leo I noted by Sohn et al. (2007) could 
be explained by collisions with Ursa Minor or Leo II, 
rather than interactions with the MW. Ursa Minor is 
also known to be kinematically disturbed (Kleyna et al. 
1998; Wilkinson et al. 2004). It has an inner bar that has 
been suggested to be tidally induced (Lokas et al. 2012) 
and has an extended stellar halo (Palma et al. 2003). To 
date, however, no sign of tidal disturbance is found for 
Leo II. 

To illustrate the non-negligible probability of such en- 
counters, we focus in Figure 14 on orbit calculations for 
the intermediate mass MW mass model (M v ; r = 1.5 x 
10 12 M©). We plot from the 10,000 randomly sampled 
points within the 4<r proper motion error space, those 
points that produce a passage of Ursa Minor (left panel) 
or Leo II (right panel) within 10 kpc of Leo I. There are 
several orbital solutions within la with such small sepa- 
rations to Leo I. Since the error space of the Leo II and 
Ursa Minor proper motions were also searched, we note 
here that solutions are also found within lcr of their re- 
spective means. So close encounters between these galax- 
ies and Leo I are not ruled out by the data, although the 
probability of such encounters is low. The encounter time 
for Ursa Minor occurs within a range of 0.8-1.1 Gyr ago, 
while the encounter time for Leo II is relatively recent, 
with a range of 0.15-0.3 Gyr ago. 



To determine whether the presence of the satellites of 
the MW can modify the orbit of Leo I, we repeated the 
analysis of Section 3.3.2, but now also accounting for the 
tidal torques exerted by the other satellites. Wc followed 
the method outlined above to test the proper motion er- 
ror space for each satellite in Monte-Carlo sense. 14 This 
showed that any multi-body tidal effects by the other 
MW satellites are insufficient to modify the average or- 
bital properties listed in Table 4. In particular, the av- 
erage velocity and time at infall are unaffected by the 
presence of the other MW satellites. We find that the 
Leo I orbit is significantly affected by the presence of 
the other satellites in only 0.1-0.2% of Monte-Carlo or- 
bits. So Leo I's high velocity is probably not a product 
of multi-body tidal torques. This makes it important to 
address whether such a high velocity can arise naturally 
in ACDM galaxy assembly scenarios without assistance 
from multi-body interactions. We explore this topic in 
detail in Paper II. 

4. CONCLUSIONS 

We have presented the first absolute proper motion 
measurement of Leo I, based on HST ACS/WFC images 
taken in two different epochs separated by ~ 5 years. 
We used the method of Sohn et al. (2012) to measure 
the average bulk motion of Leo I stars with respect to 
stationary distant galaxies in the background. We detect 
motion of Leo I at 4er confidence, and find its proper mo- 
tion to be (m^,Mjv) = (0.1140±0.0295,-0.1256±0.0293) 
mas yr _1 . The uncertainties are smaller than those ob- 
tained in previous HST studies of other MW satellites 
that used a background quasar as stationary reference. 
To derive the velocity of Leo I with respect to the MW, 

14 We decided not to explore the full proper motion error space 
of M31. The preceding analysis already established that M31 is un- 
likely to have come close enough to modify exert sufficient torques 
to modify the star formation history of Leo I. 
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we combined the proper motion with the known line-of- 
sight velocity and corrected for the solar reflex motion. 
The resulting Galactocentric radial and tangential veloc- 
ities are V rad = 167.9±2.8and V tan = 101.0i34.4kms- 1 , 
respectively. Hence, Leo I has a significant transverse ve- 
locity, but it is less than the radial velocity. 

Combined with its current position, the new knowledge 
of the three-dimensional velocity of Leo I has allowed 
us to study its orbital history in detail. To evaluate 
the past orbit we employed a range of mass models of 
increasing complexity. Starting from Keplerian models 
for the MW, we progressed first to cosmologically moti- 
vated MW models of M mr = 1.0 x 10 12 , 1.5 x 10 12 , and 
2.0 x 10 12 M Q , and then to models in which other MW 
satellites and Local Group galaxies are included as well. 
In each of these models we solved the equations of mo- 
tion to follow the Leo I velocity and position backward 
in time. We used a Monte-Carlo analysis to explore the 
impact of the observational measurement uncertainties. 

Allowing for both observational uncertainties and un- 
certainties in the MW mass model (Table 4, with a flat 
prior in mass) , we find that Leo I's most recent perigalac- 
tic passage was 1.05 ± 0.09 Gyr ago at a Galactocentric 
distance of 91±36 kpc. The ratio of the orbital pericenter 
to apocenter distance is 0.17 ± 0.07, so the orbit extends 
well outside of the MWs virial radius. Leo I entered the 
virial radius 2.33 ± 0.21 Gyr ago. This was most likely 
Leo I's first infall into the MW. A previous pericenter, 
which would have occurred almost a Hubble time ago, 
becomes slightly plausible only in the highest-mass MW 
models. 

Stellar population studies of Leo I have inferred that it 
experienced an enhancement in star formation <~ 2 Gyr 
ago. This may have been the result of Leo I entering 
the virial radius of the MW for the first time, leading 
to gas compression through ram pressure or MW tidal 
torques. Stellar population studies have also shown that 
star formation in Leo I was quenched <~ 1 Gyr ago. This 
may have been due to the pericentric approach of Leo I 
with respect to the MW, at which point ram pressure 
stripping of gas was maximal. 

A previously inferred enhanced star formation activity 
in Leo I that occurred ~ 4.5 Gyr ago is not obviously as- 
sociated with any orbital time scale, or interaction with 
any other galaxy. The orbital plane of Leo I is not coinci- 
dent spatially or in rotational sense to that of the LMC. 
The separation between M31 and Leo I remains large at 
early times, but it may have been possible that M31 ap- 
plied sufficient torques at early times to have modified 
Leo I's orbit. Leo I may have closely approached (within 
~ 10 kpc) other MW satellites, but only in the last Gyr, 
and with probabilities of at most a few percent. The 
probabilities are highest for past encounters with Ursa 
Minor or Leo II, which may have left marks in the kinc- 
matical properties of Leo I or these other galaxies. 

Given the observed velocity of Leo I and prior con- 
straints on the MW virial mass, Leo I is most likely 
bound to the MW. However, this is not true in MW 
models with masses M vir < 10 12 , at the low end of the 
allowed range. Leo I has just passed pericenter, probably 
from its first infall into the MW. So its kinematics are 
probably not virialized. So even though Leo I is probably 
bound, it is not necessarily appropriate to include it in 
equilibrium models of the MW satellite population, such 



as those that are often used to estimate the MW virial 
mass (e.g., Watkins et al. 2010). 

The velocity of Leo I can also be used to estimate the 
MW mass through the timing argument. Previous stud- 
ies of this kind have used the assumption of a radial 
orbit, but with a statistical correction for any possible 
transverse velocity (e.g., Li & White 2008). Now that 
the transverse velocity has actually been measured, it is 
possible to solve the complete timing equations for a non- 
radial orbit. This yields M M w,vir = x 1O 12 M , 
with the large uncertainty dominated by cosmic scatter. 
This is higher than, but not inconsistent with, the range 
of MW mass estimates obtained from other methods. 

This is the first paper in a series of two. In Paper II, we 
compare the new observations to the properties of Leo I 
subhalo analogs extracted from state-of-the-art cosmo- 
logical simulations. We show there that Leo I is most 
likely bound to the MW, since unbound subhalos are ex- 
tremely rare at the present epoch. We also show there 
that the observed kinematics of Leo I are more consistent 
with high-mass than with low-mass MW models. Both 
these conclusions are consistent with what we have in- 
ferred in the present paper through different arguments. 

In this paper and in Sohn et al. (2012), we have used 
our new methodologies to successfully measure with HST 
the absolute proper motions of Leo I and M31, respec- 
tively. The new measurements have allowed us to derive 
new constraints on the past and future orbital evolution 
of the target galaxies, and on the mass of the Local 
Group's dominant galaxies, the MW and M31. Moti- 
vated by these results, we are continuing to apply our 
proper motion measurement techniques to a range of 
other topics in Local Group galaxy research. For exam- 
ple, our group has ongoing HST observing programs to 
measure the proper motions of: (1) dwarf galaxies near 
the edge of the Local Group (GO-12273, PI: R. P. van der 
Marel); (2) stars at different locations in the Sagittarius 
Stream (GO-12564, PI: R. P. van der Marel); and (3) 
the dwarf galaxy Leo T, which is likely a galaxy falling 
into the MW for the first time (GO-12914, PI: T. Do). 
We expect that the results from these ongoing programs 
will further constrain the dark matter distribution in the 
Local Group and its dominant galaxies. 
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